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Regenerators  have  been  proposed  to  increase  thermal  cycle 
efficiency  of  gas  turbine  power  plants  which  are  being  considered 
for  installation  in  certain  classes  of  United  States  Navy  ships. 
These  regenerators  are  to  preheat  intake  air  and  are  currently  designed 
as  compact  heat  exchangers.   The  recommended  heat  exchanger  configur- 
ations are  either  counterflow  or  crossflow.   This  analysis  considers 
only  the  counterflow  configuration. 

Some  of  the  structural  problems  associated  with  rectangular,  com- 
pact heat  exchanger  design  are  discussed  and  certain  structural  design 
parameters  are  determined  analytically.  A  digital  computer  program  is 
presented  for  the  analysis  of  the  pressure  effects  in  an  idealized, 
two  dimensional,  rectangular  heat  exchanger.   This  program  is  written 
for  a  matrix  Analysis  of  an  idealized  heat  exchanger  with  a  cross 
section  analogous  to  a  multi-story  building  frame. 

This  investigation  was  conducted  by  Robert  L.  Corbett  and  David 
W.  Stubbs  during  the  period  from  December  1965  to  May  1966,  at  the 
United  States  Naval  Postgraduate  School,  Monterey,  California. 
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1 .    INTRODUCTION 

Several  current  projects  in  industry  and  one  at  the  United  States 
Naval  Postgraduate  School  are  actively  concerned  with  the  performance 
of  gas  turbine  regenerators.   The  impetus  for  these  investigations  stems 
from  a  proposal  to  equip  certain  types  of  United  States  Navy  ships  with 
a  power  plant  which  will  reduce,  significantly,  the  delay  time  for 
transition  from  cruising  speed  to  flank  speed. 

This  desired  flexibility  in  speed  is  planned  to  be  achieved  either 
by  augmenting  or  replacing  the  current  steam  turbine  installations  with 
diesel  or  gas  turbine  power  plants.   In  general,  the  efficiency  of  the 
diesel  is  superior  to  that  of  the  pure  gas  turbine  installations;  but 
the  advantage  in  weight  and  space  for  gas  turbine  engines  is  significant. 
These  advantages,  coupled  with  increased  efficiency  gained  through  the 
use  of  regenerators,  make  the  gas  turbine  appealing  for  shipboard 
installations.  Therefore  the  design,  operation,  and  maintenance  of 
gas  turbine  regenerators  are  important  in  the  proposed  plan  for  gas 
turbine  shipboard  installation. 

Regenerator  or  heat  exchanger  requirements  for  different  types  of 
ships  may  differ  in  size,  shape,  and  weight.  A  basic  heat  exchanger 
may  be  considered  to  consist  of  various  combinations  of  small  rectangular 
modules  stacked  or  joined  together  in  a  manner  that  allows  the  complete 
heat  exhcanger  to  meet  the  heat  transfer,  geometry,  and  weight  require- 
ments for  a  particular  installation.  The  rectangular  shape  is  a  log- 
ical choice  for  stacking  and  combining  of  modules  and  is  simple  to 
fabricate. 

Present  design  requirements  specify  maximum  pressures  up  to  ten 
atmospheres,  temperatures  ranging  from  500  to  1000  °F,  and  a  service 


life  ranging  from  20,000  to  60,000  hours.  This  latter  requirement  is 
severe,  since  structural  integrity  problems  have  occurred  in  similar 
heat  exchangers  after  operating  for  much  shorter  periods  of  time. 

The  Bureau  of  Ships  has  suggested  that  the  following  areas  be 
investigated: 

1.  Gas  side  fouling  and  air  side  salt  precipitation 

2.  Thermal  stress  behavior 

3.  Thermal  fatigue  criteria 

4.  Materials  evaluation 

5.  Header  design  and  flow  distribution 

6.  Fabrication  and  joining  techniques 

7.  Advanced  heat  exchanger  designs 

8.  Engine  evaluation  of  advanced  heat  exchangers 

This  list  indicates  some  of  the  areas  which  require  investigation  in 
order  to  define  the  requisite  design  criteria  and  to  determine  the  order 
of  reliability  to  be  specified  for  a  compact  heat  exchanger. 

Preliminary  stress  investigations  of  typical  rectangular  heat 
exchangers  have  been  accomplished  to  date.   Few  complete  analyses  of 
the  complex  state  of  stress  within  the  heat  exchanger  core  or  supporting 
frame  work  have  been  made  and  those  reported  are,  in  general,  sufficient 
only  for  the  justification  of  the  preliminary  design.   These  result 
in  heavier  and  more  costly  heat  exchanger  design,  and  possibly  in 
reduced  reliability.   The  authors  believe  that,  through  a  careful 
analytical  study  to  determine  certain  structural  design  parameters, 
heat  exchanger  design  may  be  determined  more  readily. 

The  stresses  in  different  parts  of  a  heat  exchanger  are  due 
primarily  to  the  following:  \jT\ 
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1.  The  difference  in  static  and  dynamic  fluid  pressures  and 
atmospheric  pressure. 

2.  The  thermal  expansion  or  contraction  of  component  parts 
under : 

a.  Steady  state  conditions 

b.  Transient  or  thermal  shock  conditions 

3.  Method  of  fabrication. 

The  values  of  these  stresses  are  difficult  to  determine  not  only 
because  of  the  complex  structure  but  also  since  the  value  of  pressure 
and  temperature  is  a  function  of  location  and  time. 

The  purpose  of  this  investigation  is  to  develop  a  basic  method 
to  study  the  effect  of  pressure  on  the  design.   Some  of  the  problems 
relating  to  structural  integrity  are  discussed  and  an  attempt  is  made 
to  isolate  certain  design  parameters.  The  effect  of  some  such  para- 
meters is  qualitatively  studied  through  the  use  of  the  high  speed 
digital  computer. 

The  investigation  reported  in  this  paper  was  conducted  by  the 
authors  during  the  period  from  December  1965  to  May  1966  at  the  United 
States  Naval  Postgraduate  School,  Monterey,  California.   Sincere 
appreciation  is  extended  to  the  faculty  and  staff  of  the  Aeronautical 
Engineering  Department  for  their  assistance.   The  authors  are  deeply 
indebted  to  Professor  C.  M.  Smith  for  his  invaluable  encouragement  and 
assistance  throughout  the  project.  The  authors  also  wish  to  express 
their  thanks  to  Professor  C.  H.  Kahr  for  introducing  them  to  and  in- 
structing them  in  structural  analysis  by  matrix  methods. 
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2 .    GENERAL 

Heat  exchangers  of  the  cylindrical  or  tube  type  design  have  been 
analyzed  and  have  proven  to  be  structurally  satisfactory.   On  the  other 
hand,  compact-rectangular  heat  exchangers  have  not  been  analyzed  exten- 
sively for  their  structural  action. 

A  number  of  analytical  approaches  to  the  structural  problem  were 
explored  in  this  investigation  of  pressure  effects.   Some  of  the  methods 
proved  to  be  not  entirely  satisfactory  and  sevetal  of  these  are  outlined 
in  Appendix  II. 

The  method  of  analysis  adopted  was  limited  to  the  core  and  was 
independent  of  the  flange.   The  core  is  important  from  a  heat  transfer 
point  of  view,  and  is  undoubtedly  the  most  costly.   The  flange  is  not 
included  in  this  investigation  but  its  effect  on  structural  design  is 
sufficiently  important  to  warrant  further  consideration. 

In  the  analysis  of  the  core,  cross  sections  were  considered  to  be 
similar  to  multi-story  building  frames  with  the  closure  bars  forming 
main  structural  columns,  the  fins  acting  as  partition  walls  or  columns, 
and  the  interior  plates  acting  as  floors.   Exploratory  design  parameters 
were  determined  by  considering  a  rectangular  frame  without  interior 
tie-bars  or  constraints.   The  effect  of  these  parameters  on  the  lattice- 
like frame  composing  the  core  cross  section  was  then  investigated. 
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3.   ANALYSIS 

For  the  purpose  of  this  investigation,  a  rectangular  core  module 
for  a  counter flow  heat  exchanger  was  postulated.   In  this  module  the 
geometry  and  material  properties  were  assumed  to  be  constant  over  the 
length;  the  pressure  was  assumed  to  be  constant  in  all  channels.  The 
cross  section  was  divided  into  a  pattern  of  small  rectangular  passages 
by  the  interior  plates  and  fins.  All  end  discontinuities  resulting  from 
attachment  of  the  core  to  flanges  were  neglected.   This  idealized  core 
module  is  presented  in  Figure  1. 


FIGURE  1 
IDEALIZED  COUNTERFLOW  CORE  MODULE 
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This  module  represents  only  a  small  part  of  the  complete  regenerator 
heat  exchanger.   Different  configurations  would  have  a  number  of  these 
modules  arranged  in  a  manner  suitable  to  the  operational  environment 
of  the  ship. 

Close  inspection  of  the  idealized  core  module  reveals  that  the 
fundamental  element  involved  in  the  cross  section  is  an  open  rectangular 
frame.  Many  of  these  basic  frames  are  contained  in  a  module  cross  section 
and  two  or  more  modules  are  stacked  together  to  form  a  complete  heat 
exchanger . 

Since  the  fundamental  element  is  an  open  rectangular  frame,  an 
analysis  was  performed  on  a  simple  frame  in  order  to  isolate  parameters 
which  may  affect  or  define  structural  considerations  for  design.  Qj 
In  this  analysis  the  frame  was  assumed  to  be  loaded  with  a  uniform 
internal  pressure  and  to  have  two  axes  of  symmetry  as  shown  in  Figure  2. 
While  these  assumptions  are  not  strictly  applicable  to  the  arrangement 
of  elemental  frames  in  the  core  module,  they  do  facilitate  the  analysis 
and  provide  information  as  to  which  parameters  affect  the  stresses  in 
each  member  of  the  frame. 


FZ 
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FIGURE  2 
ELEMENTAL  FRAME,  TYPICAL 
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The  two  axes  of  symmetry  which  are  shown  in  Figure  2  permit  the 
analysis  to  be  made  for  one  quarter  of  the  elemental  frame;  the  required 
boundary  conditions  are  defined  through  considerations  of  symmetry. 
The  frame  representing  one  quarter  of  the  elemental  frame  is  shown  in 
Figure  3.        L f^ 
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FIGURE  3 
ONE  QUARTER  OF  AN  ELEMENTAL  FRAME 
Considering  member  ^  of  the  frame  in  Figure  3,  the  definitive 
differential  equation  is: 


J.  A  =  a 

and  the  applicable  boundary  conditions  are: 


(1) 


a.  Ji  d3y> 
3x? 


=  o 


tf,«o 


i>.  J.d^/I  =  M, 


C.    dY,|   =  O 

d.        y,i    -  o 
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Successive  integration  of  Equation  1,  with  the  application  of 
boundary  conditions  yields: 


vcx.)    *   9*t 


M 


W 


¥ 


x?"  +  M, 


e, 


po 


6T,     J, 


(2) 


Now  considering  member  £Q   of  the  frame  in  Figure  3,  the  definitive 
differential  equation  is: 


and  the  appropriate  boundary  conditions  are : 


(3) 


dx/  k^o        2- 


b.  Oi  d2^ 


dJf^ 


x,*o 


8 


c.      dy^ 

dXi 


y1=o   ^gj,   2J; 


J.         yj      -  o 

=  0 


Integration  of  Equation  3  and  the  introduction  of  the  above 
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boundary  conditions  leads  to  the  following  expressions: 


Vfo    =    3X2.  "  Sf 


Ifyo  =  §*±  -  §££■  +  2c1-  +  m, 

v  '  2-  7.  8 


CW 


From  symmetry,  the  slope  at  the  point  X«  =  f/2  is  equal  to  zero.  This 


constraint  allows  the  determination  of  the  moment,  M^< 


M,  *  Q 


2f 


Z^-ScV-C^^) 


(5) 


-P-f  cf^) 

Assuming  that  the  sides  of  the  frame  do  not  differ  greatly  in  length 
or  in  bending  stiffness,  the  maximum  moment  occurs  at  the  corners  of  the 
frame.   From  Equation  2,  the  maximum  moment  can  be  expressed  as: 


C6) 


The  maximum  shear  force  on  the  horizontal  side  of  the  frame  is: 


Wax  -  - 


and  the  maximum  shear  force  on  the  vertical  side  is 


vm** 


=    -9£ 

"2. 


Equation  6  indicates  that  J. ,  J  ,  f/c,  and  Q  are  parameters  which 
affect  the  choice  of  design  from  a  structural  point  of  view.  Based  on 
the  assumption  that  the  pressure  and  side  dimensions  will  be  dictated 
by  requirements  for  flow  and  heat  transfer,  f/c  and  Q  may  be  fixed  for 
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a  given  design  for  reasons  other  than  structural.   Thus  the  structural 
designer  may  have  little  control  over  the  magnitude  of  the  shear  forces. 
The  magnitude  of  the  moments  may  be  controlled  by  varying  J^  and  J2. 
Therefore  the  structural  designer,  by  carefully  choosing  the  appropriate 
values  of  member  thickness,  may  be  able  to  optimize  the  design  for  the 
idealized  frame  element. 

In  order  to  study  the  effect  of  varying  member  thicknesses,  t 
and  t£,  a  two  dimensional  model  of  the  heat  exchanger  module  was  in- 
vestigated. A  diagram  of  this  model  is  shown  in  Figure  4.   This  model 
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FIGURE  4 
TWO  DIMENSIONAL  HEAT  EXCHANGER  MODEL 
shows  that  t^,  ta,  t  ,  and  tf  are  analogous  to  the  t,  and  t„  shown  in 
Figure  3.   The  effect  of  varying  t  ,  tfe ,  t  ,  and  tf  within  the  model 
chosen  is  not  obvious  and  cannot  be  defined  in  terms  of  expressions 
derived  for  moment  and  shear  for  the  simple  fr^Jne. 

In  two  dimensionalizing  the  model,  the  top  and  bottom  plates,  the 
interior  plates,  and  the  closure  bars  are  in  reality  three  dimensional 
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plates  which  have  been  Idealized  as  two  dimensional  beams  by  assuming 
an  elemental  cross  section  of  the  core.  The  Plate  Flexural  Rigidity, 
D,  is  used  for  each  element  to  include  this  effect. 

Slope  and  deflection  methods  of  analysis  are  adaptable  to  matrix 
methods  which  in  turn  are  easily  programmed  for  solution  by  the  high 
speed  digital  computer.  The  computer  program,  as  used,  is  included 
in  Appendix  I  with  a  set  of  answers  for  the  chosen  model.   This  program 
was  prepared  from  considerations  of  the  following. 

The  stiffness  method  of  matrix  solution  is  a  form  of  the  slope 
and  deflection  method  as  proposed  by  M.  L.  Pei.  [4]  The  basic  element 
for  this  method  consists  of  a  beam  under  the  action  of  end  moments  m, 
and  m2  and  end  forces  Po  and  Pa  as  shown  in  Figure  5. 
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FIGURE  5 
BASIC  ELEMENT,  SLOPE  AND  DEFLECTION  ANALYSIS 
Let  the  beam  undergo  end  rotations  0,  and  Qz  as  well  as  end  displace- 
ments S3  and  84  .  The  positive  directions  of  these  displacements  are 
as  defined  in  Figure  6. 
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FIGURE  6 
BASIC  ELEMENT  DEFLECTIONS,  SLOPE  AND  DEFLECTION  ANALYSIS 
This  beam  then,  has  four  degrees  of  freedom  represented  by  the  dis 
placement  matrix: 


W' 


Correspondingly,  loads  are  represented  by  the  load  matrix: 

four 

h  ■  te 

These  loads  and  deflections  are  related  by  the  stiffness  influence  co- 
efficient K  which  depends  upon  the  cross  section  and  the  modulus  of 
elasticity  of  the  beam.  These  equations  may  be  expressed  conveniently 
in  matrix  notation  as  follows: 


H  .  OJ  {e} 


For  the  single  member  example  given  above,  this  general  equation 
becomes : 
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^  *  KSI  ©t  f  fa^et*  K»&a  t  **4S4 


Hence  the  "K"  matrix  is  expressed  as: 


W  s 


K/i  Kiz  K13  Ki4 

Kit  Ktt  IU3  Kz* 

tot  fa  Ki&  Kb* 

k<H  Kit  K4.3  K<I4 


The  loading  IrfM     for  the  structure  is  known,  so  the  rotations  and 
deflections  of  the  member  ends,  or  more  generally,  the  joints  in  the 
structure,  may  be  determined  from: 


{6}  =  [K]"V) 


The  application  of  this  method  to  a  general  two  dimensional  frame 
consists  of  computing  the  elemental  "K"  matrix  for  each  element  of  the 
frame  and  then  adding  these  elemental  "K"  matrices  to  form  the  total  "K" 
matrix  for  the  frame.  The  "K"  matrix  is  then  inverted,  as  indicated  by 
the  above  equation,  to  determine  the  frame  joint  deflections. 

Using  the  frame  joint  deflections  the  resulting  element  loads  are 
found  from: 


{«}  -  [kJ  (©} 


where  TkJ  is  the  elemental  matrix  and  <0j  is  the  deflection  of  the 
joints  in  which  the  element  terminates.  These  elemental  loads  are  the 
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element  internal  loads.  These  internal  loads  describe  the  degree  of 
relaxation  from  that  of  a  fixed  joint.   The  external  loads  or  element 
reactions  may  then  be  determined  from: 

These  reactions  are  the  desired  results  from  which  stresses  may  be  cal- 
culated. 

The  process  of  calculating  loads  in  this  way  is  very  laborious, 
particularly  when  it  is  necessary  to  invert  a  large  "K"  matrix.   At  the 
same  time  the  process  involved  is  purely  a  mechanical  one  which  may  be 
computerized. 
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4.   RESULTS 

The  computer  was  used  to  determine  the  effect  of  varying  t  ,  t^» 

3, 

and  t  using  the  following  variations: 
P 

Basic  Configuration  Variations 

ta  =  .04  in  .04  in  --  .25  in 

tb  =  .04  in  .04  in  --  .25  in 

tp  =  .04  in  .02  in  --  .06  in 

tf  =  .01  in  Constant 

f  =  0.75  in  Constant 

c  =  0.5  in  Constant 

Maximum  normal  stresses  and  shear  stresses  for  the  variations  shown 
were  calculated.   The  shearing  stress  variations  are  functions  of  fin 
spacing  and  channel  height,  neither  of  which  was  varied  since  they  were 
assumed  to  be  dictated  by  thermodynamic  considerations.   The  shearing 
stress  variations  were  small  when  compared  to  the  normal  stress  vari- 
ations.  For  this  reason  the  variations  in  shearing  stresses  were  assumed 
to  be  negligible  for  the  parameters  being  investigated.   Graphs  showing 
the  effect  of  the  changes  on  shearing  stresses  are  presented  in  Figures 
8  and  9. 

Maximum  normal  stress  variations  for  varying  closure  bar  thickness 
are  plotted  in  Figure  10.  Maximum  stresses  are  computed  from: 

Maximum  normal  stress  -  bending  stress  plus  axial  stress. 
As  the  thickness  of  the  closure  bar  decreases; 

stress  in  the  closure  bar  is  increased  rapidly, 
stress  in  the  top  plate  is  decreased  rapidly, 
stresses  in  the  interior  plates  and  fins  are  decreased 
gradually. 
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Since  the  stresses  in  the  closure  bars  increase  rapidly,  this  method 
for  reducing  stresses  in  the  model  is  considered  marginal. 

Maximum  normal  stress  variations  for  varying  interior  plate  thick- 
ness are  shown  in  Figure  11.   This  graph  indicates  that  for  an  increas- 
ing interior  plate  thickness,  the  stresses  in  all  members  decrease 
rapidly.   However,  increasing  the  thickness  of  the  interior  plates 
rapidly  reduces  the  heat  transfer  effectiveness  between  channels.   Since 
heat  transfer  is  the  main  purpose  of  the  heat  exchanger,  this  method 
of  reducing  the  stresses  must  be  used  judiciously. 

Maximum  normal  stress  variations  for  varying  top  plate  thickness 
reduces  the  stresses  in  all  members  as  long  as  the  top  plate  is  not  too 
thick,  see  Figure  12.   Beyond  a  certain  top  plate  thickness  the  stresses 
in  the  top  plate  and  the  fins  begin  to  increase.   However,  further  thick- 
ness increase  in  the  top  plate  is  not  necessarily  detrimental  as  the  closure 
bar  stresses  continue  to  fall  off  rapidly.   At  the  point  at  which  the 
top  plate  and  the  fin  stresses  reach  a  minimum  and  the  closure  bar  stress 
reaches  a  local  minimum,  there  is  a  significant  shift  in  the  location 
of  maximum  stress  both  in  the  closure  bars  and  in  the  top  plate.   This 
maximum  stress  location  shift  results  in  a  definite  dip  in  the  maximum 
stress  values.   Since  variation  of  the  thickness  of  the  top  plate  has 
no  basic  undesirable  effects,  it  is  most  advantageous  to  make  initial 
stress  reductions  in  the  model  by  variation  of  top  plate  thickness. 

For  any  given  heat  exchanger  this  computerized  method  of  solution 
is  useful  in  the  determination  of  the  range  of  pressure  stresses  in 
the  module.   A  computer  program  and  detailed  instructions  for  its  use 
are  presented  in  Appendix  I. 
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5.    CONCLUSIONS  AND  RECOMMENDATIONS 

The  stresses  due  to  internal  pressure  are  a  significant  part  of 
the  stresses  in  a  rectangular,  compact  heat  exchanger.   The  value  of 
these  stresses  for  a  pressure  of  10  atmospheres  may  be  on  the  order 
of  sixty  to  ninety  percent  of  design  allowables. 

The  value  of  stresses  in  the  lattice-like  frame  which  composes  the 
heat  exchanger  core  can  be  controlled  through  variations  in  the  thick- 
ness of  the  individual  elements  involved. 

An  optimum  ratio  of  the  sum  of  the  stiffness  of  the  top  plate 
and  interior  plates  to  the  stiffness  of  the  closure  bars  is  implied. 
This  optimum  should  consider  heat  transfer  as  well  as  structural  re- 
quirements . 

Consideration  should  be  given  to  increasing  the  thickness  of  the 
top  plate  and/or  the  thickness  of  the  interior  plates  in  order  to 
achieve  the  required  stiffness  for  these  members.  This  approach  can 
lead  to  an  optimum  stress  level  and  distribution  for  an  internal  pressure 
loading. 

The  matrix  analysis  and  the  computer  program  used  in  this  analysis 
should  be  expanded  to  provide  solutions  for  three-dimensional  problems 
with  thermal  as  well  as  pressure  loads. 

The  compatibility  of  the  flange  and  core  is  an  important  structural 
consideration.  This  should  be  thoroughly  investigated  in  conjunction 
with  further  analysis  of  the  core. 
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FIGURE  7.  MAXIMUM  SHEAR  STRESS  VERSUS 
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FIGURE  8.  MAXIMUM  SHEAR  STRESS  VERSUS  TOP  PLATE  THICKNESS 
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APPENDIX  I 
COMPUTER  PROGRAM  "FRAME"  DESCRIPTION  AND  USE 

Program  Frame  is  a  self-contained,  Fortran  60,  digital  computer 
program  designed  to  provide  joint  deflections  and/or  rotation  and  member 
end  forces  and/or  moments  for  plane  rigid  frame  type  structures  under 
general  loading  conditions. 

The  geometry  of  the  structure  is  restricted  to  a  plane  structure 
with  all  members  parallel  to  either  the  X  or  Y  axis  in  the  usual 
cartesian  coordinate  system. 

A  generalized  flow  diagram  showing  the  basic  mechanics  of  the 
program  follows. 


Dimension  the  necessary  matrices 
AK,  f,  C,  THETA,  AKS 


Read:  NX  -  No.  degrees  kinematic  freedom. 
NM  -  No.  of  frame  members. 
EP  -  Error  parameter,  no.  author  is 
willing  to  accept  as  zero. 


Null  A  and  C  matrices. 


Read:  M  -  No.  representing  moment  deg.  of  freedom  on  left  end  of  MBR. 
N  -  No.  representing  moment  deg.  of  freedom  on  right  end  of  MBR 
K  -  No.  representing  shear  deg.  of  freedom  on  left  end  of  MBR. 
L  -  No.  representing  shear  deg.  of  freedom  on  right  end  of  MBR. 
EMOD  -  Elastic  modulus  of  the  member. 
ANER  -  Moment  of  inertia  of  the  member. 
SPAN  -  Length  of  the  member. 


© 
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Compute  elemental  AK  matrix. 


Store  AK  matrix  in  AKS  matrix, 


Store  M,  N,  K,  L  in  AKS  matrix. 


Build  A  matrix  from  AKM  matrices 


Transfer  A  matrix  into  C  matrix. 


Read:  NL  -  No.  of  loads  on  MBR. 


MBR  not  loaded 


Read:  KFM  -  Concentrated  or  distributed  load. 
W  -  Magnitude  of  external  load. 
SPA  1  -  Left  distance  from  load  to  beam  end. 
SPA  2  -  Right  distance  from  load  to  beam  end, 


Member  loaded 


1.  Replace  actual  loads  with  fixed  end  member  end  loads  and 
place  in  partitioned  "C"  matrix. 

2.  Store  fixed  end  member  end  loads  in  AKS  matrix. 


Subroutine  GAUSS  2. 

Partitioned  "C"  matrix  inverted  and  multiplied  yielding  THETA 


matrix. 


0 
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Take  elemental  "AK"  matrix  from  AKS  storage  for  one  member,  and 
multiply  times  DELTA  (THETA  for  same  member)  to  obtain  internal 
moments  and  shears . 


Add  internal  moments  and  shears  to  those  created  by  external 
loading  (stored  in  AKS  matrix)  to  obtain  external  end  moments 
and  shears. 


Print  individual  beam  moments  and  shears. 


The  following  information  is  for  those  who  wish  to  actually  use 
the  program  in  its  present  form. 

Properties  such  as  the  modules  of  elasticity  (EMOD) ,  the  area 
moment  of  inertia  (ANER) ,  and  the  span  may  vary  from  one  member  to 
another;  however,  no  provision  is  made  for  a  single  member  with  a 
variable  moment  of  inertia  (tapered  member),  except  that  an  appro- 
priate average  might  be  used  if  the  operator  so  desires. 

Introduction  of  actual  values  of  modulus  of  elasticity  and  moment 
of  inertia  will  provide  results,  in  appropriate  units,  of  actual  de- 
flections and/or  rotations.  Or,  if  desired,  the  relative  magnitude  of 
modulus  of  elasticity  and  moment  of  inertia,  based  on  selected  refer- 
ences may  be  used  instead  of  using  the  actual  values.   In  this  case, 
the  results  will  indicate  the  rotation  or  deflection  times  a  factor 
which  is  equal  to  (Eref)  x  (Iref)«   In  either  case  the  correct  values 
of  member  end  forces  and  moments  in  appropriate  units  will  result. 

The  required  input  data  is  relatively  straight  forward  and  simple 
in  that  the  first  data  card  describes  the  overall  structure  and  sub- 
sequent cards  describe  the  individual  members  taken  one  at  a  time. 
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Care  should  be  taken  in  compiling  the  data  cards,  since  seemingly 
insignificant  errors  will  produce  entirely  erroneous  results. 

The  program  provides  analysis  for  only  one  structure  at  a  time. 
This  structure  may  have  only  one  complete  set  of  loads  applied  to  it. 
The  number  of  loads  comprising  this  one  set,  however,  is  essentially 
unlimited.   Provision  is  made  for  application  of  concentrated  loads 
and/or  distributed  loads  to  any  member  and  in  any  combination. 

The  program  output  format  is  self-explanatory  except  that  THETA 
may  be  actual  rotations  or  deflections  if  actual  values  of  modulus  of 
elasticity  and  moment  of  inertia  were  used,  or  they  may  be  multiples 
of  the  reference  values  if  such  reference  values  were  used. 

Following  is  the  nomenclature  of  the  program,  a  copy  of  the  actual 
program  in  Fortran  60,  a  diagram  for  the  model  actually  solved  indicating 
member  numbers  and  degrees  of  freedom  (degrees  of  freedom  indicated 
by  arrows),  data  tabulation  for  the  model  as  shown,  a  solution  for 
deflections  and  rotations  for  the  model,  and  the  member  end  moments 
and  shears  for  the  model  and  data  cards  shown. 
Nomenclature  -  (Program  Proper) 

AK  -  Stiffness  matrix  for  an  individual  member. 
A  -  Matrix  used  for  obtaining  the  stiffness  matrix  for  the 

complete  structure. 
C  -  Stiffness  matrix  for  the  complete  structure. 
THETA  -  Matrix  of  deflections  and/or  rotations. 

AKS  -  Matrix  used  to  store  individual  AK  matrices  and  M,  N, 

K,  L,  and  end  moments  and  forces  due  to  external  loads. 
DELTA  -  THETA  matrix  for  a  single  member. 

F  -  Matrix  of  end  forces  and  moments  for  a  single  member. 
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Required  Input  Data 

Data  Card  #1   Contains  values  of  NX,  NM,  EP,  punched  in  the 
format  (2110,  E10.4). 

Data  Card  #2   Describes  the  first  member  to  be  considered. 
It  contains  values  of  M,  N,  K,  L,  EMOD,  ANER, 
SPAN,  punched  in  the  format  (412,  2E10.5, 
1F10.5). 

Data  Card  #3   Defines  the  number  of  loads  (NL)  applied  to 

the  member  described  in  card  #2.   NL  is  punched 
in  the  format  (12).   If  no  load  is  applied  to 
the  member  being  considered,  a  zero  must  be 
punched  in  column  number  2. 

Data  Card  #4   Contains  values  of  KFM,  W,  SPA  1,  SPA  2,  punched 
in  the  format  (12,  3F10.5).   If  data  card  #3 
set  NL  equal  to  zero,  this  card  (#4)  must  be 
omitted.  The  total  number  of  #4  cards  used  at 
this  point  must  equal  NL  from  card  #3. 

Subsequent  Data  Cards   The  input  data  card  sequence  is  now 
repeated  starting  with  card  #2  above.   Each 
member,  in  turn,  must  be  completely  described 
with  its  properties  and  loads  before  proceeding 
to  the  next  member. 
Nomenclature  -  (Input  Data) 

NX  -  Number  of  degrees  of  kinematic  deficiency. 

NM  -  Number  of  members  or  spans  comprising  the  structure. 

EP  -  Error  parameter,  minimum  acceptable  value  of  the  deter- 
minant of  the  stiffness  matrix.  If  determinant  is  less 
than  E2   the  program  will  stop  and  print  "matrix  singular." 
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M  -  The  number  assigned  to  the  rotational  freedom  at  the  left 

end  of  a  member. 
N  -  The  number  assigned  to  the  rotational  freedom  at  the  right 

end  of  a  member. 
K  -  The  number  assigned  to  the  translational  freedom  at  the 

left  end  of  a  member. 
L  -  The  number  assigned  to  the  translational  freedom  at  the 
right  end  of  a  member. 
EMOD  -  Modulus  of  elasticity  for  a  member.   May  also  be  used  as 

a  coefficient  of  some  reference  modulus. 
ANER  -  Area  moment  of  inertia  for  a  member.   May  also  be  used  as 

a  coefficient  of  some  reference  moment  of  inertia. 
SPAN  -  The  length  or  span  for  a  member. 

NL  -  The  number  of  external  loads  (concentrated  and/or  uniform) 
applied  to  the  single  member  being  considered. 
KFM  -  Kind  of  load  applied  (1  =  concentrated;  2  =  uniform). 
W  -  Magnitude  of  the  external  load. 
SPA  1  -  Distance  from  left  end  of  member  to  the  location  of  a 

concentrated  load. 
SPA  2  -  Distance  from  a  concentrated  load  to  the  right  end  of  a 
member . 
Basic  Element 

0  0 

ki  1l 

Note:  Vertical  members  are  rotated  90  degrees  clockwise  to  determine  the 
orientation  indicated  above. 
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The  structure  is  sketched  and  all  degrees  of  kinematic  deficiency 
are  numbered. 

The  individual  members  are  viewed  one  at  a  time,  as  shown  on 
previous  page,  and  M,  N,  K,  and  L  for  each  member  are  set  equal  to  the 
number  which  was  used  to  designate  the  corresponding  degree  of  freedom 
for  that  member. 
Load  Parameters 


KFM  =  1 


+W 


lm              SPfll               ^ 

^   spfli  _ 

1                                            <  1 

SPAM 


KFM-2 


+ 

W 
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PROGRAM  FRAME 
C  RUN  110 

READ  11»  NX.NM.  EP 
11  FORMAT  (2110.  E10.4) 

DIMENSION  AK(4»4) »A< 56.56) .C( 56.56)  »THETA(  56  )  .  AKS(  56.24  )  • 
1DELTA(4).F<4) 

DO  90  10  =  1.  NM 

READ  20.  M.  N.  K»  L.  EMOD.  ANER.  SPAN 
20  FORMAT  (412.  2E10.5.  1F10.5) 

DO  23  1-1.56 

DO  23  J=  1.56 
23  A( I »J)  =  0. 

DO  26  1*1.4 

DO  26  J«l.  4 


26 


3615 

38 

40 

42 
43 

45 
46 

4i 
49 


0 


J) 

1)*  (4.*EM0D*ANER)/SPAN 
2)*AK( 1.1) 


2) 
3) 
3) 
4) 
4) 
3) 
4) 
4) 


AK(  I 

AK(1 

AK(2 

AK(1 

AK(1 

AK(2 

AK(1 

AK(2 

AK(3 

AK(4 

AK(3 

AK(2 

AK(3 

AK(3 

AK(4 

AK(4 

AM  4 

AKS( 

AKS( 

AKS( 

AKS( 

KK=1 

DO  3615  LL=1»4 

DO  3615  JJ=1»4 

AKS( IO»KK)=AK(LL.JJ) 

KK=KK+1 

IF  (M)  103.48.38 
«  AK(l.l) 
103.  42.  40 
=  AK(1.2) 
=  AK(1»2) 
103.  45.  43 

*  AK(1.3) 
»  AK(1»3) 
103.  48.  46 
»  AM1.4) 

*  AK(1»4) 
101.  56.  49 
»  AK(2.2) 


(2«*EMOD»ANER)/SPAN 

(6.*EMOD*ANER)/(SPAN**2) 

AK(1.3) 

-AKU.3) 

-AK( 1.3) 

( 12.*EMOD*ANER)/(SPAN**3) 

AK(3.3) 

-AK(3»3) 


1)=AK( 1.2) 

1)=AK( 1.3) 

2)=AK<2.3) 

1)=AK( 1.4) 

2)=AK(2»4) 

3)=AK(3.4) 

0»17)=M 

0.18)=N 

0.19)=K 

O»20)=L 


A(M.M) 
IF  (NT 
A(M.N) 
A(N»M) 
IF  (K) 
A(M.K) 
A(K.M) 
IF  (L) 
A(M.L) 
A(L.M) 
IF  (N) 
A(N.N) 
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IF  (K)  103.53.  51 
51  A(N.K)  =  AK(2»3) 
A(K.N)  =  AK(2.3) 

53  IF  (L)  103.  56.  54 

54  A(N.L)  =  AK(2.4) 
A(L.N)  ^  AM2.4) 

56  IF  (K)  103.  61.  57 

57  A(K»K)  =  AK(3.3) 
IF  (L)  103.  61.  59 

59  A(K.L)  =  AK(3»4) 
A(L»K)  =  AK(3»4) 

61  IF  (L)  103.  63.  62 

62  A(L»L)  =  AK(4>4) 

63  DO  65  1=1.56 
DO  65  J=1.56 

65  C( I.J)  -    C( I.J)  +  A( I.  J) 

NPM=NX  +1 

READ  681.  NL 
681  FORMAT  (12) 

IF  <NL)  103.  90.  69 
69  DO  90   KO  =  l.NL 

READ  71.  KFM.  W»  SPA1,  SPA2 
71  FORMAT   (112.  3F10.5) 

IF  (KFM-2)  73.  82.  103 

73  IF  (M)103.741.74 

74  C(M.NPM)  =  (W*SPA1*SPA2**2)/(SPAN**2)  +  C(M.NPM) 
741  AKS( I0»21)=-(W»SPA1*SPA2**2)/(SPAN**2)  +AKS(I0.21) 

IF  (N) 103.761.76 
76  C(N.NPM)  =  -(W*SPA2*SPA1**2)/(SPAN**2)  +  C(N.NPM) 
761  AKS( I0»22)=(W*SPA2*SPA1**2)/(SPAN**2)  +  AKSU0.22) 
IF  (K)  103.781.78 
78  C(K.NPM)  =  ( (W*SPA2**2)*(3.*SPA1  +  SPA2 )  )  / (SPAN**3 )  +  C(K.NPM) 
781  AKS( IO»23)=-( ( W*SPA2**2 )* ( 3.*SPA1+SPA2 ) )/(SPAN**3)  +AKS( 10.23) 
IF  (L)  103.801.80 
80  C(L.NPM)  =  ( (W*SPA1**2)*(3.*SPA2  +  SPA1 ) ) / ( SPAN**3 )  +  C(L.NPM) 
801  AKS( I0»24)=-( (W*SPA1**2)*(3.*SPA2+SPA1) )/(SPAN**3)  +AKS( 10.24) 
GO  TO  90 

82  IF  (M)  103.831.83 

83  C(M.NPM)    =    (W*SPAN**2)/(12»)    +    C(M.NPM) 
831    AKS( I0.21)=-(W*SPAN**2)/(12. )    +AKSU0.21) 

IF    (N)     103.851.85 
85    C(N.NPM)    =    -(W*SPAN**2)/(12. )    +    C(N.NPM) 
851    AKS( I0.22)=(W*SPAN**2)/(12.)    +AKS(I0.22) 
IF(K>103»871,87 
87    C(K.NPM)    =    (W*SPAN)/(2«)    +    C(K.NPM) 
871    AKS(IO»23)=-(W*SPAN)/(2.)    +AKSU0.23) 
IF(L)     103.891.89 

89  C(L.NPM)    =    (W*SPAN)/(2« ) ■+    C(L»NPM) 
891    AKS(I0#24)=-(W*SPAN)/(2»)    +AKS(I0»24) 

90  CONTINUE 
NP  *  1 

£Abb  6AU992  ?NX»  NP»  EP#  G»  THETA.  Kl  ) 
GO  TO  (93.  99)  .  Kl 
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93  PRINT  94 

94  F0RMAT(////44H1N0TE   1.  THETA  IS  JOINT  ROTATION/DEFLECTION/) 
PRINT  321 

321  FORMAT  (29H  NOTE   2.  ANSWER  *  !I  X  THETA////) 
DO  96  K=1»NX 

96  PRINT  97,  K»  THETA(K) 

97  FORMAT  (11H0EI  X  THETA  I3»3H  =  F20.15) 
PRINT  9072 

9072  F0RMATC////43H1MEMBER  END  MOMENTS  AND  SHEARS  APPEAR  BELOW/) 

PRINT  9074 
9074  FORMAT(47H  CLOCKWISE  MOMENTS  AND  DOWN  SHEARS  ARE  POSITIVE////) 

PRINT  730 

730  FORMAT( ////43H  LEFT      RIGHT      LEFT      RIGHT) 
PRINT  731 

731  FORMAT(43H  MEMBER    MOMENT    MOMENT     SHEAR     SHEAR//) 
DO  9735  IM=1»NM 

KJ=  1 

DO  976  LL=1,4 

DO  976  JJ=lt4 

AK(LL,JJ)=AKS( IM»KJ) 
976  KJ=KJ  +1 

M=AKS( IM»17) 

IF  (M)  103.9711,979 
979  DELTAt 1)=THETA(M) 

GO  TO  9712 

9711  DELTA(1)=0. 

9712  N=AKS( IM,18) 

IF(N)  103, 9716,9714 
9714  DELTA(2)=  THETA(N) 
GO  TO  9717 

9716  DELTA(2)=0. 

9717  K=AKS( IM»19) 

IF  (K)  103,9721,9719 
9719  DELTA(3)  =THETA(K) 
GO  TO  9722 

9721  DELTA(3)=0. 

9722  L=AKS( IM,20) 

IF  (L)103,9726,9724 
9724  DELTA(4)=THETA(L) 
GO  TO  9727 

9726  DELTA(4)=0. 

9727  DO  9730  KB=1,4 
F(KB)=0. 

DO  9730  JB=1»4 
9730  F(KB)=F(KB)  +AK(KB,JB)  *  DELTA(JB) 
F(1)=F(1)  +AKS(IM,21) 
F(2)=F(2)  +  AKS( IM,22) 
F(3)=F(3)  +   AKS(IM,23) 
F(4)=F(4)  +  AKS( IM»24) 

9735  PRINT  9736t  IM,F( 1 ) , F ( 2 ) »F ( 3 ) ,F ( 4 ) 

9736  FORMAT  (I4»  F9#3»3F7,3) 
60  TO  105 

99  PRINT  100 
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100 

103 
104 
105 


10 


11 
12 
13 


14 
20 
30 
31 


32 

33 
34 
40 


41 

42 

43 
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FORMAT  (1 
GO  TO  105 
PRINT  104 
FORMAT  (1 
CONTINUE 
STOP 
END 

SUBROUTI 
DIMENSION 
NPM=N+M 
DO  34  L* 
KP  =  0 
Z  =  0.0 
DO  12  K= 
IF(Z-ABS 
Z=ABSF(A 
KP  =  K 

CONTINUE 
IF(L-KP) 
DO  14  J= 
Z=A(L»J) 
A(L,J)=A 
A(KP.J)= 
IF(ABSF( 
IF(L-N)3 
LP1=L+1 
DO  34  K= 
IF(A(K.L 
RATIO=A( 
DO  33  J= 
A(K.J)=A 
CONTINUE 
DO  43  1= 
II=N+1-I 
DO  43  J= 
JPN=J+N 
S  =  0.0 
IF(II-N) 
IIP1-II+ 
DO  42  i* 
S  =  S+A(  II 
X(  II)  =  ( 
KER  =  1 
RETURN 
KER  =  2 
END 
END 


6H  MATRIX  SINGULAR) 


2H  INPUT  ERROR) 


NE         GAUSS2(N.M»EP.A.X.KER) 
A(56»56) »X(56) 

1»N 


LtN 

F(A(K.L) ) Ul.12.12 

(K.L) ) 


13.20,20 
LfNPM 

(KP.J) 

Z 

A(L.L) )-EP)50»50»30 

1»40»40 

LPltN 

)  )32»34.32 

K.L)/A(L.L) 

LP1.NPM 

(K.J)-RATIO*A(L,J) 

1»N 

1»M 


41.43.43 

1 

IIP1.N 

»K)*X(K) 

A( II»JPN)-S)/AUI»II ) 
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'* 

>* 
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V 
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^ 

* 
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N 
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N 

<K 

** 

>a 

,v 

V 

> 

'V 
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^ 

ft. 

* 

<K 

'J* 

N 

t 

? 

§ 

y 

y 

V 

V 

* 

* 
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'* 

V 
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t 
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1020039 

1 

2      -] 

2033940 

1 

2      -] 

3044041 

1 

2      -] 

4054100 

1 

2      -] 

6070042 

>2 

51. 2E 

55    .: 

01     .IE 

LE-07 
01 

.75 

.50. 
51. 2E 

01 

0. 
•  IE 

0. 
01 

.75 

.50. 
51. 2E 

01 

0. 
•  IE 

0. 
01 

.75 

.50. 
51. 2E 

01 

b. 

.IE 

0. 
01 

.75 

.50. 

•  8E 

01 

0. 
.IE 

0. 
01 

.375 

7084239 

.8E 

01 

.IE 

01 

.375 

8093943 

.8E 

01 

•  IE 

01 

.375 

9104340 

•  8E 

01 

.IE 

01 

.375 

10114044 

•  8E 

01 

.IE 

01 

.375 

11124441 

•  8E 

01 

.IE 

01 

.375 

12134145 

•  8E 

01 

•  IE 

01 

.375 

13144500 

•  8E 

01 

•  IE 

01 

.375 

15160042 

•  8E 

01 

•  IE 

01 

.375 

16174246 

•  8E 

01 

.IE 

01 

.375 

17184643 

•  8E 

01 

.IE 

01 

.375 

18194347 

•  8E 

01 

•  IE 

01 

.375 

19204744 

•  8E 

01 

.IE 

01 

.375 

20214448 

.8E 

01 

•  IE 

01 

.375 

21224845 

•  8E 

01 

•  IE 

01 

.375 

22234500 

•  8E 

01 

•  IE 

01 

.375 

24250049 

•  8E 

01 

•  IE 

01 

•  375 

25264946 

•  8E 

01 

.IE 

01 

•  375 

26274650 

•  8E 

61 

•  IE 

ei 

.375 
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27285047  .8E  01 

28294751  .8E  01 

29305148  .8E  01 

30314852  .8E  01 

31325200  .8E  01 

33340049   51. 2E  01 

1 

2  150. 
34354950   51. 2E  01 

1 

2  150. 
35365051   51. 2E  01 

1 

2  150. 
36375152   51. 2E  01 

1 

2  150. 
37385200   51. 2E  01 

1 

2       150. 

6010000     .8E  01 

1 

2      -150. 

8020000  .0125E  01 

10030000  .0125E  01 
12040000  .0125E  01 

14050000     .8E  01 

1 

2  150. 
15060000     .8E  01 

1 

2  -150. 
16070000  .0125E  01 

18090000  .0125E  01 

20110000  .0125E  01 

22130000  .0125E  01 


.IE  01      .37"5 

•IE  01      .375 

•IE  01      .375 

.IE  01      .375 

•IE  01      .375 

•IE  01      .375 

0.  0. 

•IE  01       .75 

0.  0. 

•IE  01       .75 

0.  0. 

•IE  01       .75 

0.  0. 

•IE  01       .375 

0.  0. 

.IE  01       .5 

0.  0. 

.IE  01       .5 

•IE  01       .5 

•IE  01       .5 

•IE  01       .5 

0.  0. 

.IE  01       .5 

0.  0. 

.IE  01       .5 

•IE  01       .5 

•IE  01       .5 

•IE  01       .5 
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23140000 

1 
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24150000 
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2 
26170000 
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01 

.IE 

01 

.5 
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•  IE 

01 

0. 

•  5 
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.5 

28190000 

.0125E 
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•  IE 

01 

.5 

30210000 

•0125E 

01 

.IE 

01 

.5 

32230000 

1 

2 
33240000 
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2 
34250000 

•  8E 

01 

•  IE 
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.5 
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•  8E 

01 
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•  IE 
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0. 

.5 
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.0125E 
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i 

•  IE 

01 

0. 

.5 

35270000 

.0125E 

01 

.IE 

01 

•  5 

36290000 

•0125E 

01 

•  IE 
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•  5 

37310000 

•0125E 

01 
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01 

.5 
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01 

•  IE 
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150. 
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-ftm-E fr—THETA    IS— JfM-NT-RO  T  AT  10  N  /DEFLECT  I  ON 

NOTE   2.  ANSWER  =  EI  X  THETA 


EI  X  THETA 

'E'l— X— THETA- 

EI  X  THETA 


THETA 
THETA   5  = 

THE-T* 6— =" 

THETA  7  = 
THETA  8  = 
THETA  9  = 
THETA— 10-=- 
THETA  11  = 
THETA  12  = 
THETA  13  = 
THET*— YH-=- 
THETA  15  = 
THETA  16  = 
THETA  17  = 


18  = 

19  = 

THETA  20  = 
THETA  21  = 

-£-I— X— TH  E  T  A-^2-^- 
EI  X  THETA  23  = 
EI  X  THETA  24  = 
EI  X  THETA  25  = 

"Ei— X— TH  E  T  A-2-6— =- 
EI  X  THETA  27  = 
EI  X  THETA  28  = 
EI  X  THETA  29  = 

-E  I— X— T  H  E  T  A— 3  0-f 
EI  X  THETA  31  = 
EI  X  THETA  32  = 
EI  X  THETA  33  = 

-El— X— TH  E  T  A— 5  U— =" 
EI  X  THETA  35  = 
EI  X  THETA  36  = 
EI  X  THETA  37  = 

-E  I— X— TH  E  T-A-  3  8~^- 
EI  X  THETA  39  = 
X  THETA  40  = 
X  THETA  41 


EI 
EI 

-El 
EI 
EI 
EI 

-El- 
EI 
EI 
EI 

-El- 
El 
EI 


-THET-A-U-2-^ 
THETA  43  = 
THETA  44  = 
THETA  45  = 
-THE-T-A— 4-6— =" 
THETA  47  = 
THETA  48  = 
THETA  49  = 
-TH  E-T-A-50-^- 
THETA  51  = 
THETA  52  = 


.059757 
.039902 
.000000 
.039902 
.059757 
V020913 
.084237 
.000181 
.030079 
-.000000 
.030079 
.000181 
.084237 
.020913 
.023057 
.038969 
.041627 
.012380 
.000000 
.012380 
.041627 
.0  38969 
.023057 
.025138 
.003051 
.053146 
.01  1316 
.000000 
.01 1316 
.053146 
.003051 
.025138 
.049U25 
.047713 
.018272 
.018272 
.047713 
.04  9425 
.039716 
.054848 
.039716 
.005222 
.019661 
.019661 
.005222 
-.01  1389 
.012368 
.01 1389 
.018564 
~rO  4-4 199 
.044199 
.018564 


1681  41914 

351722048" 

000003752 

351729324 

1681 16448 

817393648 

065282650 

155194702 

360703894 

00000682  1- 

360741638 

155228802 

065302659 

817406391- 

530810092 

665726654 

451475536 

074253315- 

000023057 

074174139 

451528257 

665821242 

530988808 

163619431 

132390510 

538107285- 

232965555 

000938599 

233752268 

539722547 

131717825 

164648979 

681714638 

020395349- 

389301728 

389324466 

020370793 

681679168- 

653672258 

966682534 

653668620 

375323569- 

487614485 

487617668 

3753247D6 

984892958- 

798751 140 

984899324 

758140656 

61496296-3" 

614952049 

758128377 
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M E  H  B  C R    E  Nfl-MflMEW~S— AMD-  SH€A  R  S~ A  P P  E  A  R-~8 E  LO  W~ 


CLOCKWISE    MOMENTS    AND    DOWN    SHEARS    *RE    POSITIVE 


MEMBER 


LEFT 
MOMENT 


MOMENT 


LEFT 
SHEAR 


RIGHT 
SHEAR 


1 
2 

5 


9 
ID 
11 
-1-2- 
13 
14 
15 
4  6- 
17 
18 
19 


-20- 
21 
22 
25 

-2^4- 


25 
26 
27 
-2-8- 
29 
30 
31 


-3-2- 
53 
34 
35 


-36- 
37 
38 
39 

-40- 
Ul 
42 
43 

-4V 
45 
46 
47 

-4.6- 


S-T-QP- 

TIMi, 


49 
50 
51 
-5-2- 
53 
54 
55 


6.28 

19.29 

•21.13 


-+9.-33- 

-.03 
4.58 
-8.14 
-9-.-44- 


•10.73 

9.42 

-8.  16 

— 4-.-5Jh 

1.48 

-4.12 

7.57 

— 9-.-S8- 

10.40 

-9.88 

7.56 

—4,-1-2- 

-4.06 

4.98 

-7.15 

-9-.-90- 


•10.38 

9.90 

-7.  11 

—5.00- 

-7.09 

8.20 

17.92 


+7-.-93- 


-7 


25 

.55 

.02 

.00- 

.02 

.55 

2 .  32 

—.00- 

-.00 

.00 

.00 

—2r.-32- 


4.00 
.03 

-.00 
.-0-3- 


•4.00 

7.09 

.05 

-r0+ 


-.01 

-.05 

■7.09 


19 
21 
19 
—6 
-4 
8 
-9 


33 

13 

29 

28- 

51 

16 

4? 


+0V73- 


-9 

8 

-4 


4 

-7 

9 


U4 

14 

58 

.-03 


12 

56 
88 


— T0-.40 
9.88 
-7.57 
4.12 

-K48- 

-5.00 

7.11 

-9.90 

— 10-.-38- 

-9.90 

7.15 

-U.98 

4.06- 

-8.25 

-17.93 

-17.92 

— -8,20 

7.09 

-6.28 

-.04 

-*00 


.04 
6.28 
-2.52 

.06- 

-.02 

.02 

.06 

— 2-.-52- 

-3.80 

-.02 

-.00 

— .02- 

3.80 
.07 
.03 

— nrGO- 


.00 
-.03 
-.07 


90.39 

58.71 

53.79 

—22.11- 

-12.1  1 

33.98 

-46.84 

^-3-.-7^- 


-53.79 
46.84 
-33.98 
— 1  2.-11- 
14.94 
-31.15 
46.54 
—54.08- 
54.08 
-46.54 
31.15 
—  14.94- 
-24.17 
32.24 
-45.45 
—54.0  8- 
-54.08 
45.45 
-32.24 
—24.1-7- 
-69.05 
-69.22 
-56.25 
—43.26- 
12.80 
30.04 
-.12 

.00- 

.12 
-30.04 
37.09 
—,14- 


-.05 

.05 

.14 

-3T-.09- 

37.90 

.03 

-.00 

,03- 


-37.90 

51.82 

.15 

,02- 

-.02 

-.15 

-51.82 


22.11 

53.79 

58.71 

-90.39- 

12.11 

•33.98 

46.84 

■53.79 

53.79 

•46.84 

33.98 

■12,11 

•14.94 

31.15 

•46.54 

-54.08- 

•54.08 

46.54 

•31.15 

-14.94- 

24.  17 

•32.24 

45.45 

-54.08- 

54.08 

•45.45 

32.24 

-24.17- 

12.80 
■43.28 
•56.25 
69.22- 
•69.05 
44.96 
.12 
— .00 
-.12 
•44.96 
37.91 
,  1  4- 


.05 

-.05 

-.  14 

-37-.T01- 

37.10 

-.03 

.00 

.-03- 

•37.10 
23.18 
-.15 
— ^r02- 


.02 
.15 
23.18 
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APPENDIX  II 
OTHER  METHODS  OF  ANALYSIS 

Prior  to  the  decision  to  use  the  two-dimensional  frame  analysis 
for  the  internal  pressure  stresses  a  number  of  other  methods  were  in- 
vestigated.  Generally,  the  approach  was  one  of  searching  for  the  correct 
design  philosophy  and  attempting  to  adapt  simple  solutions  to  the  prob- 
lem at  hand. 

The  methods  outlined  herein  seemed  to  apply  to  the  problem  but  in 
each  case  the  results  became  too  complex  to  provide  useful  design  in- 
formation or  the  results  required  specific  knowledge  of  some  particular 
parameter  which  was  unavailable. 

TOP  PLATE  ON  AN  ELASTIC  FOUNDATION 

The  first  approach  was  to  consider  the  top  plate  of  the  heat  ex- 
changer to  be  fixed  on  the  two  sides  adjoining  the  side  closure  bars, 
to  be  resting  on  an  elastic  foundation,  with  modulus  K  lb/in,  along  its 
span,  and  to  be  loaded  with  a  uniform  pressure  load,  Q  lb/in.   The  elastic 
foundation  was  assumed  to  be  the  result  of  the  arrangement  of  interior 
plates  and  fins.   Assuming  no  change  in  geometry  or  other  parameters 
with  length,  the  plate  was  analyzed  as  a  beam  considering  the  stiffness 
to  be  EI/(1-1>2)  instead  of  simply  EI. 

The  expression  for  the  moment  at  the  end  of  this  beam  is: 


2A2  LSNH^US.ryXbJ 


where:  AHFTZ^n 


**    V  4£I 


This  result  shows  that  the  modulus  of  the  elastic  foundation,  K,  is  of 
primary  importance.   In  truth,  this  modulus  is  a  function  of  such  para- 
meters as  fin  spacing,  width  of  the  heat  exchanger,  and  thickness  or 
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stiffness  of  the  interior  plates.  No  adequate  analytical  means  of 
estimating  the  magnitude  of  the  foundation  modulus  was  found.   It  is 
recommended  that  K  be  determined  by  testing  actual  modules. 
DISCONTINUITIES  BETWEEN  HEADER  AND  CORE 
A  number  of  attempts  were  made  to  determine  the  value  of  discon- 
tinuity stresses  at  the  joint  between  the  header  or  flange  and  core. 
In  general,  this  method  attempted  to  extend  the  classical  solutions  for 
cylindrical  pressure  vessels  to  include  rectangular  pressure  vessels. 
This  approach  might  be  possible  if  the  corners  of  the  rectangular  heat 
exchanger  could  be  rounded  off  slightly  such  that  a  continuous  mathematical 
function  could  be  used  to  define  the  surface  of  the  heat  exchanger.   The 
best  approach  with  the  square  corners  appears  to  be  with  the  use  of  an 
elemental  lengthwise  strip  taken  at  the  location  of  the  greatest  dis- 
continuity in  deformation  between  the  core  and  the  flange  as  shown  below. 

Ft.flNfiE 


If  QQ  and  Mq  could  be  determined,  the  elemental  strip  of  the  top 
plate  of  the  core  could  be  analyzed  as  a  beam  on  an  elastic  foundation 
with  end  loads  QQ  and  M  .   Difficulties,  in  connection  with  finding  QQ 
and  Mq,  the  magnitude  of  the  elastic  foundation  modulus,  whether  the 
elemental  strip  of  the  core  could  be  considered  semi-infinite  in  length, 
and  the  correct  method  for  determining  the  d  between  the  header  and 
core,  were  encountered. 
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FRAME  ANALYSIS 

A  number  of  lattice-type  frames  were  analyzed  and  the  analysis  for 
the  simple  frame  with  two  axis  of  symmetry  did,  in  fact,  provide  the 
impetus  for  continuing  with  the  complex  frame  analysis  as  used.   This 
analysis  provided  fundamental  information  as  to  the  effect  of  the  various 
parameters  involved  in  the  core. 

Frames  with  various  numbers  of  non-yielding  tie  bars,  which  repre- 
sented the  interior  plates,  were  analyzed.   The  solutions  neglected  the 
effect  of  fins.   This  approach  was  abandoned  because  of  the  inability 
to  adequately  account  for  the  effect  of  the  fins. 

SANDWICH  CONSTRUCTION 

Considerable  effort  was  expended  in  researching  publications  dealing 
with  sandwich  construction.   Much  is  available  on  the  subject  in  the 
literature  but  no  reference  was  found  which  dealt  with  sandwich  or  honey- 
comb structures  loaded  with  an  internal  pressure. 


50 


INITIAL  DISTRIBUTION  LIST 


No.  Copies 


1.  Defense  Documentation  Center  20 

Cameron  Station 
Alexandria,  Virginia  22314 

2 .  Library  2 

U.  S.  Naval  Postgraduate  School,  Monterey,  Calif. 

3.  Chairman,  Department  of  Aeronautics  2 

U.  S.  Naval  Postgraduate  School,  Monterey,  Calif. 

4.  BUSHIPS  (Code  645)  1 

Department  of  the  Navy 
Washington,  D.  C.   20360 

5.  Prof  C.  M.  Smith  1 

Aeronautical  Engineering  Department 

U.  S.  Naval  Postgraduate  School,  Monterey,  Calif. 

6.  LT  R.  L.  Corbett  1 

HELANTISUBRON  THREE 

FPO,  New  York,  New  York  09501 

7.  LT  D.  W.  Stubbs  1 

1224  Big  Horn  Ave. 
Alliance,  Nebraska 

8.  Prof  C.  H.  Kahr  1 

Aeronautical  Engineering  Department 

U.  S.  Naval  Postgraduate  School,  Monterey,  Calif. 

9.  Commander,  Naval  Air  Systems  Command  1 

Department  of  the  Navy 
Washington,  D.  C.   20360 


51 


UNCLASSIFIED 


Security  Classification 


DOCUMENT  CONTROL  DATA  •  R&D 

(Security  classification  of  title,  body  of  abstract  and  indexing  annotation  must  be  entered  when  the  overall  report  la  etaaallled) 


I.  ORIGINATING  ACTIVITY  (Corporate  author) 

U.    S.    Naval    Postgraduate  School 


2a.    REPORT  SECURITY    CLASSIFICATION 

UNCLASSIFIED 


26     CROUP 


3.  REPORT  TITLE 

An  Investigation  of  the  Pressure  Stresses  in  a  Compact,  Plate  Fin, 
Counter  flow  Heat  Exchanger 


4     DESCRIPTIVE  NOTES  (Type  ot  report  and  Inclusive  datea) 

Thesis 


5.  AUTHORfSJ  (Last  name,  first  name,  initial) 

Corbett,  Robert  L. ,  LT,  USN 
Stubbs,  David  W. ,  LT,  USN 


6-  REPORT  DATE 

May  1966 


7a.     TOTAL   NO.    OF     PAGES 


50 


7  b.    NO.    OF   REFS 


8a      CONTRACT    OR    GRANT    NO. 


6.    PROJECT  NO. 


9a.    ORIGINATOR'S   REPORT   NUMBERfS) 


9b.  OTHER  REPORT   NOfS)  (Any  other  numbera  that  may  be  assisted 
(his  report) 


This  document  has  been  approved  for  public 
release  and  sale;  its  distribution  is  unlimited. 


10.  AVAILABILITY/LIMITATION  NOTICES 


11.  SUPPLEMENTARY  NOTES 


12.  SPONSORING  MILITARY  ACTIVITY 

Commander,  Naval  Air  Systems  Command 
Department  of  the  Navy 
Washington,  D.  C,  20360 


13.  ABSTRACT 

Regenerators  have  been  proposed  to  increase  thermal  cycle  efficiency 
of  gas  turbine  power  plants  which  are  being  considered  for  installation  in 
certain  classes  of  United  States  Navy  ships.   These  regenerators  are  to 
preheat  intake  air  and  are  currently  designed  as  compact  heat  exchangers. 
The  recommended  heat  exchanger  configurations  are  either  counter  flow  or  '■ 
cross  flow.   This  analysis  considers  only  the  counter  flow  configuration. 

Some  of  the  structural  problems  associated  with  rectangular,  compact 
heat  exchanger  design  are  discussed  and  certain  structural  design  parameters 
are  determined  analytically.   A  digital  computer  program  is  presented  for 
the  analysis  of  the  pressure  effects  in  an  idealized,  two  dimensional, 
rectangular  heat  exchanger.   This  program  is  written  for  a  matrix  analysis 
of  an  idealized  heat  exchanger  with  a  cross  section  analogous  to  a  multi- 
story building  frame. 

This  investigation  was  conducted  by  Robert  L.  Corbett  and  David  W. 
Stubbs  during  the  period  from  December  1965  to  May  1966,  at  the  United 
States  Naval  Postgraduate  School,  Monterey,  California. 


nn  form 

\J  \J       1  JAN  64 


1473 


UNCLASSIFIED 


53 


Security  Classification 


UNCLASSIFIED 


Security  Classification 


14- 


KEY  WORDS 


LINK  A 


HOLE 


LINK  B 


ROLE 


WT 


LINK  C 


Heat  Exchanger 
Matrix 
Rectangular 
Stress 


INSTRUCTIONS 


1.    ORIGINATING  ACTIVITY:    Enter  the  name  and  address 
of  the  contractor,  subcontractor,  grantee,  Department  of  De- 
fense activity  or  other  organization  (corporate  author)  issuing 
the  report. 

2a.    REPORT  SECURITY  CLASSIFICATION:    Enter  the  over- 
all security  classification  of  the  report.    Indicate  whether 
"Restricted  Data"  is  included.    Marking  is  to  be  in  accord- 
ance with  appropriate  security  regulations. 

2b.    GROUP:    Automatic  downgrading  is  specified  in  DoD  Di- 
rective 5200. 10  and  Armed  Forces  Industrial  Manual.    Enter 
the  group  number.    Also,  when  applicable,  show  that  optional 
markings  have  been  used  for  Group  3  and  Group  4  as  author- 
ized. 

3.  REPORT  TITLE:     Enter  the  complete  report  title  in  all 
capital  letters.    Titles  in  all  cases  should  be  unclassified. 
If  a  meaningful  title  cannot  be  selected  without  classifica- 
tion, show  title  classification  in  all  capitals  in  parenthesis 
immediately  following  the  title. 

4.  DESCRIPTIVE  NOTES:    If  appropriate,  enter  the  type  of 
report,  e.g.,  interim,  progress,  summary,  annual,  or  final. 
Give  the  inclusive  dates  when  a  specific  reporting  period  is 
covered. 

5.  AUTHOR(S):    Enter  the  name(s)  of  author^)  as  shown  on 
or  in  the  report.    Enter  last  name,  first  name,  middle  initial. 
If  military,  show  rank  and  branch  of  service.    The  name  of 
the  principal  author  is  an  absolute  minimum  requirement. 

6.  REPORT  DATE:    Enter  the  date  of  the  report  as  day, 
month,  year,  or  month,  year.    If  more  than  one  date  appears 
on  the  report,  use  date  of  publication. 

7  a.    TOTAL  NUMBER  OF  PAGES:    The  total  page  count 
should  follow  normal  pagination  procedures,  i.e.,  enter  the 
number  of  pages  containing  information. 

7b.    NUMBER  OF  REFERENCES:     Enter  the  total  number  of 
references  cited  in  the  report. 

8a.    CONTRACT  OR  GRANT  NUMBER:    If  appropriate,  enter 
the  applicable  number  of  the  contract  or  grant  under  which 
the  report  was  written. 

8b,  8c,  flt  8d.  'PROJECT  NUMBER:  Enter  the  appropriate 
military  department  identification,  such  as  project  number, 
subproject  number,  system  numbers,  task  number,  etc. 

9a.    ORIGINATOR'S  REPORT  NUMBER(S):    Enter  the  offi- 
cial report  number  by  which  the  document  will  be  identified 
and  controlled  by  the  originating  activity.    This  number  must 
be  unique  to  this  report. 

9b.  OTHER  REPORT  NUMBER(S):  If  the  report  has  been 
assigned  any  other  report  numbers  (either  by  the  originator 
or  by  the  sponsor),  also  enter  this  number(s). 

10.    AVAILABILITY/LIMITATION  NOTICES:    Enter  any  lim- 
itations on  further  dissemination  of  the  report,  other  then  those 


imposed  by  security  classification,  using  standard  statements 
such  as: 

(1)  "Qualified  requesters  may  obtain  copies  of  this 
report  from  DDC  " 

(2)  "Foreign  announcement  and  dissemination  of  this 
report  by  DDC  is  not  authorized. " 

(3)  "U.  S.  Government  agencies  may  obtain  copies  of 
this  report  directly  from  DDC.    Other  qualified  DDC 
users  shall  request  through 


(4)      "U.  S.  military  agencies  may  obtain  copies  of  this 
report  directly  from  DDC    Other  qualified  users 
shall  request  through 


(5)     "All  distribution  of  this  report  is  controlled.    Qual- 
ified DDC  users  shall  request  through 


If  the  report  has  been  furnished  to  the  Office  of  Technical 
Services,  Department  of  Commerce,  for  sale  to  the  public,  indi- 
cate this  fact  and  enter  the  price,  if  known. 

1L    SUPPLEMENTARY  NOTES: 
tory  notes. 

1Z  SPONSORING  MILITARY  ACTIVITY:  Enter  the  name  of 
the  departmental  project  office  or  laboratory  sponsoring  (pay- 
ing for)  the  research  and  development.    Include  address. 

13.  ABSTRACT:    Enter  an  abstract  giving  a  brief  and  factual 
summary  of  the  document  indicative  of  the  report,  even  though 
it  may  also  appear  elsewhere  in  the  body  of  the  technical  re- 
port.   If  additional  space  is  required,  a  continuation  sheet  shall 
be  attached. 

It  is  highly  desirable  that  the  abstract  of  classified  reports 
be  unclassified.    Each  paragraph  of  the  abstract  shall  end  with 
an  indication  of  the  military  security  classification  of  the  in- 
formation in  the  paragraph,  represented  as  (TS),  (S),  (C),  or  (V). 

There  is  no  limitation  on  the  length  of  the  abstract.    How- 
ever, the  suggested  length  is  from  150  to  225  words. 

14.  KEY  WORDS:    Key  words  are  technically  meaningful  terms 
or  short  phrases  that  characterize  a  report  and  may  be  used  as 
index  entries  for  cataloging  the  report.    Key  words  must  be 
selected  so  that  no  security  classification  is  required.    Identi- 
fiers, such  as  equipment  model  designation,  trade  name,  military 
project  code  name,  geographic  location,  may  be  used  as  key 
words  but  will  be  followed  by  an  indication  of  technical  con- 
text.   The  assignment  of  links,  roles,  and  weights  is  optional. 


Use  for  additional  explana- 


DD 


FORM 

1    JAN  04 


1473  (BACK) 


UNCLASSIFIED 


54 


Security  Classification 


?l7^TfrilllVjfvffi^^ilf^ffl'jfl^CTlfMH, 


thesC7546 

An  investigation  of  the  pressure  stresse 


3  2768  002  09151  4 
DUDLEY  KNOX  LIBRARY 


InHHi 


WWmMm 


SSSfittJ 

$S! 

lllilll 

■;,;;).l:;: 

He! 

»> 

uSSHAi 


.■;:f.>3^ffl» 


